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Abstract. Given an increasing function H : [0, 1) —> [0, oo) and 



A n (H) := w£ \ J2j t ' (*< - t)H 2 (t)d?j 



where T„ := {r = (ti)f =0 : = to < t\ < • ■ ■ < t n = 1}, we characterize the property 
A n {H) < and give conditions for A n (H) < and A n (H) > -^=p for /3 € (0, 1), 
both in terms of integrability properties of H. These results are applied to the approxi- 
mation of certain stochastic integrals. 
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I. Introduction 

In this paper we estimate the size of the error which occurs when a 
stochastic integral is approximated discretely. To explain the problem 
in more detail, we assume a stochastic process X = (X t )te[o,i] such that 

dX t = a(X t )dW t with X = x > 0, 

where W = {W t )te[Q,i} is the standard Brownian motion, a satisfies 
certain regularity properties, and (J 7 t)te[o,l\ is the augmentation of the 
filtration generated by W. It is of interest to approximate discretely a 
stochastic integral, which can be written as 

(1) f(X 1 )=Ef(X 1 )+ [ \ u dX u , 

Jo 

where / : M. — > R is a polynomially bounded, Borel measurable func- 
tion, and A = (A t ) t e[o 1] is a suitable adapted process. We approximate 

/(■V;) l.v 

n 

(2) E/(^) + 5]Vi^-U 

1=1 

where r n := (tj)™ =0 is a deterministic time net with = t < t\ < 
• • • < t n _i < t n = 1. Using this approximation instead of the original 
stochastic integral, we obtain an approximation error 

rt 

(3) fix,) - E/(X0 - K-A*u - Xt^). 

1=1 
1 
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We are interested in the minimal quadratic error under the constraint 
that the time net used in the approximation has n + 1 time points. 
According to [SI Lemma 3.2 and its proof], this error is equivalent to 

(4) <£(Z) :=inf a x (Z,r n ), 

where a x (Z, r n ) := inf (E|/(X 1 ) - Ef(X 1 ) - ^( X U ~ ' 
and Z = f(Xi) with the infimum taken over all sequences v = (fj)^T 
of jF^-measurable step functions Vi : Q — > R. 



The approximation problem is of interest for at least two reasons. 

(a) In stochastic finance one would like to replace a continuously 
adjusted hedging portfolio in the Black-Scholes option pricing 
model by a discretely adjusted one, as portfolios can be adjusted 
in practice only finitely many times. If we consider the qua- 
dratic error which occurs in this replacement (and which we can 
interpret as risk in finance), then we end up with the approxi- 
mation problem described above. In this case X = (X t ) tg [ 01 ] is 
an appropriate positive diffusion process, / : (0, oo) — > [0, oo) is 
a payoff function of a European type option, and r n is the net 
of time points where the portfolio is rebalanced. 

(b) The approximation introduced above yields to an approxima- 
tion of j'l \ u dX u by YJi=\ K-i ( x um - X u _ lM ). The point is 
that the approximation itself is a stochastic integral, but the 
integrand X u (which is usually hard to compute) is only com- 
puted n-times, whereas the increments (X tiM — X^m) can be 
easily simulated (for example by an Euler scheme). 

There are several previous results concerning the error caused by the 
discrete approximation of stochastic integrals. Under certain conditions 
on Z and a, C. and S. Geiss showed that if r n = (^)" =0 is the equidistant 
time net with cardinality n + 1, then one has that 

ax(Z,r n )< ' 
'n 



if and only if Z belongs to the Malliavin Sobolev space Di^ [3[ The- 
orems 2.3, and 2.6]. Furthermore, they proved that there exists a 
constant c > such that a*(Z) > -^-^ unless there are constants c 
and ci such that Z = Co + c\X\ a.s. [3> Theorem 2.5] (if such constants 
do exist, then a*(Z) = 0). It is also known by (3J Theorem 2.9] and [3, 
Theorem 3.2] that there exists a constant c > such that a*(Z) < 
if Z has a certain polynomial smoothness measured by Besov spaces 
generated by real interpolation. In this case the rate is obtained by 
adapted non-equidistant time nets. 
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M. Hujo showed in [HI Theorem 3] that for X being the Brownian 
motion or the geometric Brownian motion, there exists random vari- 
ables Z = f(Xt) G L 2 (0,.F,P) such that 

sup y/na^(Z) = oo, 

which means that the approximation rate is not always ^= even if the 
underlying process is the standard Brownian motion. However, there 
are no explicit examples of such functions. 

These results lead us to the question to characterize those Z = 
f(X 1 ) G L 2 (tt,F,F) with 

a n(Z) — ~~i= f° r some c = c(Z) > 0. 



According to Theorem 14.41 below, there exists a constant c = c(a) > 
such that 

1 / n rU \ 5 

-a x (Z, r)<[22J (k - t)H 2 x Z{t)dt J < ca x (Z, r), 



J=l "H-l 



where H x Z{t) := || (<x 2 0) ^ X t)\ L , F : [0,1) x I ^ R is given by 

F(t, x) = K(Z\X t = x), f and X satisfy certain conditions and Jcl 
depends on X. Moreover, Lemma [4.31 implies that H X Z is increasing 
so that we concentrate our investigations for some time on the quantity 

A n (H) := mf [J2 [ ' (U~ t)H\t)dt 

" \i = l J^- 1 / 

where the function H : [0, 1) — > [0, oo) is increasing and 

% : = {r = (ti)? =0 : = t < *i < • • • < t n = 1}. 



Our first main result, Theorem 12. 3[ says that 



inf (j2 [ {ti-t)H 2 (t)dt 
reTn \i=i Jti-i , 



i 

2 

< C 



n 



if and only if the function H is integrable. Moreover, in Theorem 12.51 
we give sufficient conditions for 



i 

2 

c 



and 

i 



H (E /' (ti-t)H 2 (t)dt) >^= 
where (3 G (0, 1), in terms of the growth rate of H : [0, 1) — > [0, oo). 
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These results can be applied to the setting introduced above and 
also to other situations, for example to the quadratic approximation of 
multi-dimensional stochastic integrals (see [9], [12] and [TB"]). 

The paper is organized as follows: In section 2 we introduce the main 
results of the paper, their proofs can be found in Section 3. In Section 
4 we apply the results of Section 2 to the 1-dimensional stochastic 
setting. In particular, we give an example of random variables for which 
the approximation rate is -^=^ < a*{Z) < ^=y, for (3 G (0, 1) in case 



X is the standard Brownian motion or the geometric Brownian motion. 
In Section 5 the results of Section 2 are applied to the approximation 
of d-dimensional stochastic integrals where the underlying diffusion 
might have a drift. 

2. Results 

To shorten the notation in the following, we say that A ~ c B if for 
constant c > 1 it holds that \A < B < cA and for time net r 6 T n we 
define 

| |t| |oo := max {U - U-x}. 

i£l,...,n 

2.1. Definition. Let H : [0, 1) — > H. be a non-negative measurable 
function. If r = (tj)™ =0 G T n , then we define 



A n {H) :=inf r6rn A(F,r). 

2.2. Definition. We say that an increasing function H : [0, 1) — > [0, oo) 
belongs to the set A if and only if 

:= sup y/nA n (H) < oo, 
and to the set 7i if and only if 

\\H\\ n := [ H{t)dt < oo. 



2.3. Theorem. Let H : [0, 1) — > [0, oo) be an increasing function. Then 

sup y/nA n (H) < oo 



if and only if f H(t)dt < oo. In particular, one has that 

\\H\\a \\H\\h- 

2.4. Remark. The proof of Theorem 12.31 implies that / := H{t)dt < 
oo gives 

A n (H) < -1= for all n e N. 



n 
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This rate can be obtained by regular sequences (see [TU] and [TT] ) 
generated by H. Regular sequences generated by H are time nets 
T n = (£")f =0 for which 

r f ? % r 1 

/ H(t)dt = - / H(t)dt 
Jo n J 

for all i G {0, n}. 

Our second main result is 

2.5. Theorem. Let H : [0,1) — > [0, oo) fre an increasing function and 
a G (|, 1). T/ien one /ias i/ie following: 

(1) J/ £/iere exists a constant c% > 1 sitc/i i/ia£ 

(1 - loefl - 01 _a 
H(t)< Cl y Sl _ J; /or aZ/t G [0,1), 

i/ien 

A„(if) < - ° /or a// fiGN, 



2a- 1 



n 

where c = c(a) > 1. 
(2) J/ t/jere exists s G [0, 1) and a constant c 2 > 1 such that 

g(f)s l (i-i°g(i-t))- /oroii(eM)i 

c 2 (i — tj 

AJif) > - ; 1 /or all nGN, 
where c = c(s, a, Cq) > 1. 



2.6. Remark. It follows from the arguments in [SJ Lemma 4.14, Propo- 
sition 4.16] that if H is increasing and there are C G (0, oo), a G (1, oo) 
with 

. * 

(t) " [a + log(l + ^)]-(l-t) 
for all t G [0, 1), then one has that 

sup y/nA n (H) < oo for all nGN. 

n 

2.7. Remark. Let H : [0,1) — ■> [0, oo) be a measurable function and 

:= sup sg r o t i H(s) < oo for all t G [0, 1). Then the monotonicity 
properties of A n (-) imply the following: 

(1) sup ngN ^JnA n {H) < \\H*\\ n as a consequence of Lemma 3.1. 

(2) if H*{t) < c ^f-vr" for all t e [ 0; 1); then 



c 



AJH) < - for all nGN. 



n 
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3. Proof 

In this chapter we prove Theorems 12.31 and 12.51 To prove Theorem 
12 .31 we need two lemmas concerning the connection between A n (H) and 
/„ H{t)dt, where H is a non-negative and increasing; function. 

3.1. Lemma. Let H : [0, T) — > [0, oo), T > 0, be an increasing function 
such that 



1=1 H(t)dt < oo. 
Jo 



Then for all n G N there exists a sequence r" = (£™)™ =0 ; = t^ < t\ < 



• • • < t™ = T such that 



1 * H(t)dt = -I 
Jo n 

for all i < n and for this sequence it holds that 



Proof. The existence of the sequence (£™)™ =0 for which 



1 H(t)dt = —I 
n 



follows from the continuity of the integral. Now we have 

/ ' (t?-t)H\ t )dt = j2 / * m-t)H(t)]H(t)dt 

1=1 Jt i-1 i = l J h-i 

<-J2 sup ($-t)H(t). 

n i=l *£[*?-!.*?) 

Since H is increasing, it is clear that 

(%-t)H(t)< / H(s)ds< / H(s)ds 
Jt Jt^, 

for all t G [*?_!,*"). Hence 

V / (U- t)H 2 (t)dt <~J2 H{t)dt = —. 
i=i ^-i n i=i ^-i n 

□ 

3.2. Lemma. Lei : [0,1) — > [0, oo) fre an increasing function. If for 
all n G N t/iere exists a tzme net r n = (t™)™ =0 G 7^ such that 

A(H, r") < 4= 
\/n 
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for some fixed c > 0, then H is integrable and 



[ H(t)dt < V2c. 
Jo 



Proof. If A(H, T n ) = 0, then H = and the claim is trivial. Assume 
then that A(H, r n ) > 0, which implies that H(t) > for some t G [0, 1). 
Let a := inf{t G [0, 1) : H(t) > 0} and f n = {a} U {tj» G r n : t? > a}. 
Since H is positive on (a, 1), our assumption implies that Ht™!^ — > 
as n — > oo. Using the Cauchy-Schwartz inequality and the assumption 
A 2 (#, r n ) < ^ we see that 

2 



(5) 



n-l 



n-l 



i=l 

n . t n 

< 2n V / (t? - t)H 2 {t)dt < 2c 2 . 
i=i ^-i 

Let b G (a, 1) and < e < a/c. Choose n such that b < and 

„fo n-l 

/ ff(t)dt<^ff(^ 1 )(t?-Ci) + e. 

i=l 

(We can choose n satisfying this, since the positivity of the function H 
on the interval (a, 1) implies that t^„ 1 — ► 1 and | \f n \ \ 00 — ► as n — > oo.) 
Now (jSJ) gives that 

/ H(t)dt < V2c + e 
Jo 

and since this is true for any b G (a, I) and any e > 0, we finally have 

H(t)dt < V2c. 



f 

Jo 



□ 

Proof of Theorem 12. 31 Assume first that H E Ti.. Then / := H{t)dt < 
oo and Lemma [3. II implies 

y/nA n (H) < I for all n G N 

and < ||-^||w- 

Assume now that if G .A, which means that 

sup y/nA n (H) < oo. 

n&i 

Lemma 13.21 implies 

H(t)dt < v^sup ^nA n {H) 



n&i 

and \\H\\ n < y/2\\H\\ A . 
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The computations above imply 



/ H(t)dt < oo if and only if sup y/nA n (H) < oo 
Jo «eN 

and that ll-^IU- D 

3.3. Lemma. Let H : [0, 1) — > [0, oo) 6e an increasing function. Then 

i 

2 



A n (#) < inf 
V y ~ Te(o,i) 



(lo T H(t)dt) rl 

± -^-+ / (1 - t)i^ 2 (t)^ 

n — 1 



/or all n > 2. 



Proof. Let T e [0, 1) and let r n = (tj)" =0 G T n be a time net such that 
= t < ti < ■ ■ ■ < t n -i = T < t n = 1 and 

"U a rT 



H(t)dt 



n — 1 



H(t)dt for all i — 1, n — 1. 



'o " x 

Using Lemma [3. II we get that 

A 2 (H,r n ) = J2 (U-t)H 2 (t)dt+ I (l-t)H 2 (t)dt 



i J ti-i 



tn-l 



< 



I T H(t)dt 



n 



-^-+ / (l-£)# 2 (t)dt. 



By definition, we have that A n (H) < A(H, r n ) and we are done. □ 

3.4. Remark. The best rate that Lemma 13.31 can give, is obtained by 
choosing T such that 



J H(t)dt = v^^T^y {l-t)H 2 (t)dt 



1/2 



However, it is not known if Lemma 13.31 gives the optimal rate, i.e. we 
do not know whether the inequality 

2 



(6) A 2 n (H)>- inf 
c re(o,i) 

holds. What we have is 
A 2 n (H)=Jnf 

where 



Jo H ^) dt 



n 



at rl 

-^-+ / {l-t)H 2 (t)dt 
1 Jt 



Te(0,i) 



A 2 n _ l (H\[0,T]) + £(l-t)H 2 (t)dt 



n-l ft 



4-i(#l[0,T]) := inf V / ( tl -t)H 2 (t)dt. 

U=tO<— <*n-l=J . — T ./( 
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In order to obtain inequality © we would need to know that there 
exists a constant c > such that 

, 2 i 2 



AUM0,T])> c n _ 1 , 

for all n > 2, but we do not know whether this is true. 

For the proof of Theorem I2.5[ we need the following lemmas. 

3.5. Lemma. Let (3 G (0, 1). Then there exists a constant c > such 
that 

(1 - log(l - t))-( 1+ V foo 



2 



/oo 
z-P- 2 (l-t)$- 2 dz for all t G [0,1). 



Proof. Let ^(t) = (IzMlzglZ^ and ^( t ) = /» ^(1 _ t) \^ dz _ 
Choosing x = — Iog ^~^ ; we obtain 

/oo 
z -f>-\\-t) k .- 2 dz 

/ dog(l-t)\^ 2 /i . . 2 log(l-t)da; 

— — (1 - t) — — it- 1 — 

X J X 1 



log(l-t) 

-log(l-t))-^ 1 r-^-t) 



x^e x dx, 



since (1 - 1) ^s(i-t) = [ e i°g(i-*)]io g (i-t) = e . 
The statement follows from 



lim ^44 = I" x p e~ x dx G (0, oo). 



□ 



3.6. Lemma. 0, Lemma 7] Let 9 G [1,2) and if e : [0, 1) — > [0, oo), &e 
given &?/ 

H e {t) = ^(2-9)(l-t)- 6 forte [0,1). 

TTien 

n . 
/ H 

inf 

(*■ 

/or a// n G {1, 2, ...} 



inf V / - > (0 ~ l) n_1 



3.7. Lemma. Lei Lf : [0,1) — > [0, oo) be an increasing function and 

/3g(o,i). if 

/oo 
z- - 2 (l-t)^ 2 dz for all te [0,1), 

£/ien 

A„(Lf) > — -= for allneN 

cpvn 13 
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where c p = ^ f3(¥+ 2 + 2?+ 2 + l)e. 
Proof. Let g : [1, oo) x [0, 1) — > (0, oo) be given by 

g(z,t) = z-e~\l-t)^. 

Then 

^ (M) = (l + -Y +2 (i-t)^ <2^+ 2 
g(k + l,t) \ + k) { } 

for all k > 1 and i G [0, 1). We have 

^(s, t) = (- log(l - t) - (2 + P)z) {l ~Z~ 2 

and it is easy to see that for any fixed t G [0, 1) there exists k t > 2 such 
that g(z, t) is increasing for all z < k t — 1 and decreasing for all z > k t . 
Hence 

g(z,t)dz> ^g(h,t)+ 9(k,t), 

k=l k=k t +l 

where we treat an empty sum as zero. Since g(k,t) < 2 l3+2 g{k + l,t) 
for all k > 1, we have 

0(A* - 1, f) + + ^(fct + 1, t) < c p g(k t + l,t), 

with C/3 := (4 /3+2 + 2 /3+2 + 1), and therefore 

oo oo 

9(k,t) = g(k t -l,t) + g(k t ,t) + g(k t + l,t) + 9(k + l,t) 

fc=(fc t -l) k=k t +l 

oo 

<c p ^2g(k + l,t). 

k=kt 

This implies 



poo k t~2 oo 
1 f i 7 



&=1 k=kt 
kt — 2 ^ oo 



>^^(A;,*) + - <?(M) 
fc=i 0/3 fc=(fc t -i) 

>-Ej(m) 



for all * G [0,1). 
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Let Qk = 2 — j and pk = k~( 1+/3 \ By assumption, 



H 2 (t) > j g(z,t)dz 

- oo 

-5>(M) 



oo 

1 

> 

~~ C /3 



fc=l 
oo 



S- 2 



1 A 1 1 _ l 

p k=i 
1 v 

-)] ft (2-a t )(l-{p. 



Now 



A n( H ) =i £j2 {U-t)H\t)dt 

i=i 

^itM («.-«)-5>e- 

1=1 P fc=l 

P fc=l 8=1 • /t »-l 



a fc )(l-t)- afc ^ 
-a k ){l-ty a *dt 



> — y^p* h / ft - ^)( 2 - a fc)( i - *r afe ^. 



p k=l i=l Jt *-i 



To prove our claim it is enough to consider n > 2. We set 
H ak (t) = y/(2-a k )(l-t)-^ 



and now Lemma [3.61 implies that 

1 oo 

AH(fl)>-^K-ir j 



oo X 1 X n_1 



-(1+0) 



C /3 t~f V *: 



fc=l 
1 oo 

^ fc=n 

1 



1 



where cp = eficp. □ 



12 



HEIKKI SEPPALA 



3.8. Lemma. Let (3 G (0,1) and H : [0,1) — > [0, oo) be an increasing 
function such that there exists a constant c\ > 1 /or which 



A n (H+l)> 



civn 



for all n£N. 



Then there exists a constant c 2 > 1 such that 

1 



A n (if) > 



C 2 Vfl 13 



/or a// n e N. 



Proof. Assume first n > n := {2 l3+1 c\) 1 ~f ) . Then we have 2 c L {2ny — n 
and since 



Al n _ 1 (H + l)<2[Al(H)+Al(l)]<2 
we get 



for all n e N, 



> 



> 



2c 2 (2n-l)/ 3 2n _ 4c 2 (2n)^ cfn/ 3 



for all n > n, where c 2 = 2 t \[2c\. 

If n < h, the computations above imply 



where C2 = C2 |~n] 2 and [n] := inf {k 6Z:ii< A;}. 



□ 



Proof of Theorem 12.51 

(1) Let T = l-e c « (n) , where c a (n) = 1 - ((1 - a)n l - a + l)i==. Then 



< ci 



(l-log(l-t))- Q 



1-t 



dt 



ci 



1 — a 

cm 1 "" 



[(1 — log(l — T)) — 1] 



and 



(1 - t)H 2 (t)dt < c\ 



^l-Ml-t))- 2 " 
1 - t 



-dt 



2a - 1 

r 2 

2a - 1 

-.2 



L-(l-log(l-T)) 



l-2a 



((l-cyK^ + l)- 



< C l(l - a) 1 -° l-2q 

2a- 1 
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and hence Lemma [3.31 says that, for n > 2, 
1 



A n (H) < 



n — 1 



T 



H(t)dt) 



+ 



< 



n — 1 



1 2- 2a i 2~ l-2ct 

n + qc a n 



1 

T 

1/2 



[l-t)H\t)dt 



1/2 



< c r 



f2 + c Q )2 



n 



2a- 1 



l-2a 

(l-a)~r=5T 
2a-l 



where c c 

(2) Assume there exists a constant C2 > 1 such that 



H(t) > (1 " l0 f (1 ", t)r " foralUG[ S ,l). 



c 2 (l - 1) 

Then there exists a constant C3 > 1 such that 

;i-io g (i-t))- Q 



H{t) + 1 > 



for all t e [0, 1) 



ca(i-t) 

If we write /3 = 2a — 1 G (0,1), Lemma [3.51 implies that there exists a 
constant C4 > 1 such that 

1 /"°° 

{H{t) + 1) 2 >- z- p - 2 (l-t)-*- 2 dz for alltG [0,1), 
C4 Ji 

and Lemma 13.71 implies that there exists C5 > 1 such that 

1 



4.0ff + i) > 



for all nGN. 



C5V nf 3 

Finally, Lemma T3.8I implies the existence of a constant c > 1 such that 

1 



A n (H) > 



for all nGN. 



cv n h 



□ 



4. Application: Optimal approximation rate of certain 

stochastic integrals 

Throughout the section, we assume a standard Brownian motion 
W = (W t )t<E[o,i] on a stochastic basis (ft, P, (Tt) te [o,i]), where (J^) te [o,i] 
is the augmentation of the natural filtration of W and T = T\. We 
let the process S = (S , i ) t6 [ 0i i] be the geometric Brownian motion, i.e. 
St = e Wt ~^ for all t G [0, 1]. The space of continuous, infinitely many 
times continuously differentiable functions with bounded derivatives is 
denoted by C£°(R). Moreover, we let X = (Xt)te[o,i] be a diffusion such 
that 



(7) 



dX t = a(X t )dW t with X = x G R, 
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where the process X is obtained through Y = (Vf)t e [o,i] given as unique 
continuous solution of 

dY t = a{Y t )dW t + b(Y t )dt with Y = y e R, 

with < e < b G (R) and 6 G C£°(R), in the following two ways: 

(a) y = £o e R, a := a, b := 0, X t := Y t , 

(b) y = logx with x > 0, 

a(y):=^l, b(y) := -\a{y)\ and X t = e Y \ 

Moreover, we let 7 be the Gaussian measure on R, i.e. 

1 _^ 
d r y(x) := e 2 c?x. 

V 27T 

4.1. Definition. Let C e be the linear space of Borel measurable func- 
tions / : R — > R such that there exists m > for which 

supe- mN E/ 2 (x + t#) < 00 

for all t > 0, where (7 is a centered standard normal random variable. 
Moreover, we define 

C := {Z := f(Y{) :il^R\feC e and Y as above}. 

The main tool for investigating the approximation problem in papers 
of C. Geiss, S. Geiss, and Hujo was the iJ-functional defined in the 
following way. 

4.2. Definition. Let X be a stochastic process as in (J7J) and assume 
that Z G C (or Z G L 2 (Q, T, P) if X G {W, 5}). Then we set 



(8) iJ x Z(t) := 



d 2 F 



for all t G [0, 1), 

L 2 



where F : [0, 1) x J -> R is given by = E(Z|X t = x), with 

/ = R in the case of (a) and / = (0, 00) in the case of (b). 

4.3. Lemma. [3j Lemma 5.3], [3, Lemma 3.9] The function HxZ : 
[0, 1) — > [0, 00) is continuous and increasing. 



In order to deduce from Theorem 12.31 a characterization of the ap- 
proximation rate 



we need the following theorem. 
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4.4. Theorem. [3J Lemma 3.2] [7J Lemma 3.10] Let X be a stochastic 
process as in Z G C (or Z G L 2 (Q,.F,P) if X G and 
r=(t i )^ =0 GT n . Then 



i=l 



a x (Z,r) ~c(J2 {U-t)H 2 x Z{t)dt 



where c > 1 is an absolute constant depending on a only. Consequently, 

A n (H x Z). 

D and Z G C for Z G L 2 (0,.F,P) i/ 



4.5. Corollary. Let X be as in 

X G {W,S}). Then 



sup Vna^ (Z) 

nGN 



(7 



,d 2 F 
dx 2 



(t,x t 



dt. 



L 2 



with I 



where F : [0, l)xJ->lis given by F(t, x) = E(Z\X t = x) 
in the case of (a) and I = (0, oo) in the case of (b). 

Proof. Theorem 12.31 together with Lemma 14.31 and Theorem 14.41 gives 
the result immediately. □ 



2 d 2 F 
dx 2 



{t,x t 



is inte- 



4.6. Remark. Remark 12.41 implies that if 

grable, then the regular sequences generated by (^°' 2 ^^j 
give the rate -7=. Using these sequences, denoted by r™, we have that 



dt < 00, then 



< gig for all nEN, 



n 



a*{Z)<a x {Zy r 



where cg^j > is taken from Theorem 14.41 above. 

One can also optimize over random time nets instead of deterministic 
ones considered here. The result [U Theorem 1.1.] from C. and S. Geiss 
implies that ^= is the best possible approximation rate also for the 
random time nets in case the underlying diffusion X is the Brownian 
motion W or the geometric Brownian motion S and Z is not equal to 
Co + C\X\ a.s. for some Co, c\ G K. This means that if X G {W, S}, the 
random time nets do not improve the approximation if the deterministic 
time nets already give the rate -^=. According to this, Corollary 14.51 
implies that if 



a 



,d 2 F 
dx 2 



(t,x t 



dt < 00, 



then the optimal approximation rate is also for the random time 
nets and this rate is obtained by using the regular sequences generated 



by 



a 



2 d 2 F 
dx 2 



(t,x t 
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Now we give for (3 G (0, 1) an example such that 

a%(Z) ~ c -= for all fiGff, 

in case X is a standard Brownian motion or the geometric Brownian 
motion. According to Theorem 12.51 Lemma 14.31 and Theorem 14.41 it is 
sufficient to find a random variable Z = f a (Wi) such that 

Exm „. (1-^(1-'))- , 

where a = ^p. 

4.7. Example. Let a G (1/2,1) and f a = Y^=o a khk G £2(7), where 

a = (a k )T= is s iven b y 

3 if k G {0,1,3}, 

0, = S 75 [ik 2 > 

k-2 



and (/ifc)feLo ^- -^2(7) is the complete orthonormal system of Hermite 
polynomials, 

h k (x) = ' e 2 —re 2 . 
Then Z a := G L 2 (Q, J 7 , P) and it can be shown that 

1 /2 

ff^«(t) = ( 1 + ^ifctog-^ ) ~ C1 (1 ~ 1 ° 1 g( ^" t)ra 



fc=2 



for all t G [0, 1) (according to Lemmas 14.91 and 14.81 below). Using 
Lemma 14.91 it is easy to show that there exists a constant C2 > such 
that 

H w Z a (t) ~ C2 H s Z a (t) for all t G (0, 1). 
Theorem 12.51 implies there exists a constant C3 > 1 such that 

i <at(Z a )<-^= 



for all n G N, where X G {W, S}. In other words, letting (3 G (0,1) 
2 



and defining a := pp we have 



a*(Z n ) — = for all n G N. 



The following lemma should be known. For completeness and con- 
venience of the reader we include a proof. 
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4.8. Lemma. Let f3 > 1. Then for all t G [0, 1), one has that 

where the constant c > 1 depends at most on (5. 

Proof. Let n > e 13 be an integer, e G [^-,^), and t = e~ e . Since 
k\og~ f3 (k) is increasing if k > e 13 and we assumed that n > e 13 , we have 

oo 2n 

1 + J2 kl °Z~^( k ) tk >^k\og~P{k){e- l/n ) k 

k=2 k=n 
In 

k=n 

> e~ 2 n 2 log _/3 (n). 

Moreover, 

oo n oo (m+l)n 

1 + J]A;log- /3 (A;)t fc < 1 + ^ fclog^(fc) + ^ ^ /c log~ /3 (/c)e~^ 

fc=2 fc=2 m=l k=mn+l 

n oo 

< cg^nlog _/3 (n) + + l)n 2 log _/3 (n)e - ^ 

fc=2 m=l 

oo 

< C/3 n 2 log-^(n) + n 2 log _/J (ra) ^(m + l)e" m/2 

m=l 

< (c/3 + c)ra 2 log -/3 (n), 

where depends at most on /3 and c = Ylm=i( m + l)e~ m//2 . This 
implies, for t = e~ e with e G [^j, that 

oo 

1 + J2 k ^g~ p (k)t k ~ C1 n 2 log _/J (n) for all n > e^, 

fe=2 

where ci > 1 is a constant depending at most on (5. Adapting the 
constant C\ > 0, we get this for n > 2. 
Now we show that if n > 4, then 

(i ~yy^ ~- »'■■**(")■ 

where c 2 > 2 is a constant depending at most on /3. Firstly, we have 
that log(j) ~ C3 ^, where c 3 = |. Moreover 

log(w _1 ) ~ C4 1 it, 

for all if G [e -1 / 2 , 1], where c 4 = [2(1 — e~5)] _1 . Hence 

1 

1 -t ~ CB -, 

n 
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where C5 = |[1 — e^^]^ 1 . Furthermore, 

lncr ri 

< log(n/c 5 ) < log((l — t) ) < log(c 5 n) < 21ogn 



2 

since C5 < 2 and n > 4. Now 

1 - log(l - t) ~ 3 log(n 

and hence 



(1 - log(l - t))-f 2 _ p 

~„- n log M n, 



where c 2 = 3^0^ . 



(1-t) 

5 



C2 



If t > e 3 , the computations above imply that 

(1 ~ lT (1 ~ t)r ' ~ C2 - 2 log* n ~ Cl 1 + £ Jfe log-<W, 

where n is such that e » < t < e n + l . If < t < e 4 , then one has 
that 



1 < 1 + J2 kl °S~ l3 (k)t k < Cp, 



k=2 

where the constant cp > depends only on /3, and 
1 (l-log(l-t))-^ 

^- - d/ * 

where the constant dp > depends only on /3. Hence 

^ ' fc=2 

for all t G [0, 1), where the constant c > 1 depends on /3. □ 

4.9. Lemma. [7J Lemma 3.9] For / = Y^k=o a khk € ^2(7), £ € [0, 1) 
and Z = /(Wi) one /ias t/iat 



^Z 2 (t) = ^a2 +2 (fc + 2)(fc+l)t fc , 



fc=0 

= E ( a *+ 2 - ^ff J (* + 2)(* + l)t*, 

where W is a standard Brownian motion and S is the geometric Brow- 
nian motion. Moreover 

^H w Z 2 (t) - jj(a 2 + a 2 ) < H s Z\t) < AH w Z 2 {t) + 2a 2 . 
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5. Application: Approximation of certain d-dimensional 
stochastic integrals with drift 

We can apply Theorems 12.31 and 12.51 also to the discrete time ap- 
proximation of (i-dimensional stochastic integrals considered by Zhang 
[13] . Temam [12] and Hujo [9]. Our setting introduced below recalls 
for the convenience of the reader line by line the setting of [9], which 
generalizes the 1-dimensional setting of Section 4 to d dimensions. 

We assume a stochastic basis (f2,.F, P, (J-t)te[o,i]), where (J r t )te[o,i] is 
the augmentation of the natural filtration generated by the ci-dimensional 
standard Brownian motion W = (W t ) t e[o,i] with T = T\. 

We consider a diffusion X = (X 1 , ...,X d ), where 

nt d ft 

(10) X\ = 4+ h{X u )du + V / <Tij{X u )dWi t G [0, 1], a.s. 

Jo j=1 Jo 

for all i = l,...,d and Xo = (xq,...,Xq). We assume X is obtained 
through Y given as unique path-wise continuous solution of 

ft d ft 

(11) F/ = ^+/ b t (Y u )du + J2 Vij( Y u)dWi, te [0,1], a.s. 

Jo j=1 Jo 

for alH = 1, d, where b h G C£°(R d ) and {aa T ) ij (x) = J2i=i &ik(x)&jk 
is uniformly elliptic, i.e. 

d 

(v°- T )ij(x)&j > A||CH 2 , for a11 x,^eR d and some A > 0, 

where || • || is the Euclidean norm. Again, we assume that X is obtained 
through Y by one of the following two ways: 

(a) x = y e R d , bi(x) := bi(x), <7y(x) := <Tij(x), and X t = Y t , 

(b) x = G (0,oo) d , k{y) := ^? - §E-=i<^)> *M : = 

Here and in the following e y = (e Vl , ...,e Vd ) for y = (y%, y^). As in 
one dimensional case, (a) is related to the standard Brownian motion 
and (b) to the geometric Brownian motion. 

Moreover, we assume that / : E — > R is a Borel-function such that 
for some q G (0, oo) and C > it holds that 



(12) \f(x)\ < C(l + \\x\\ q ), xeE 

where the set E is defined by 

l d in case (a) 

(0, oo) d in case (6). 



£ := 
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Finally, we define the function g : M. d — > R by 

f(y) in case (a) 



g(y) ■-- 



f(e y ) in case (b). 



5.1. Theorem. (TJ Theorem 8 on p. 263], [2, Theorem 5.4 on p. 149] 
For b, a with aa T uniformly elliptic, there exists a transition density 
T : (0, l]xl rf xR d ^ [0, oo) G C°° such that F(Y t G B) = f B T(t, y, £)d£ 
fort G (0, 1] and B G B(M. d ), where Y = (F t ) tg [ 0) i] is the strong solution 
of stochastic differential equation f TTTT) starting in y. Moreover, the 
following is satisfied: 

(i) For (s,y,£) G (0, 1] x R d x R d we have that 

d 1 d d d 2 d d 

k,l=l j=l yK yt i=l yt 

(ii) For a G {0,1,2,...} and multi-indices b and c there exist positive 
constants C and D, depending only on a, b, c and d such that 

Qa+\b\ + \c\ 



< C „- D M^Ml 



t {d+2a+\b\+\c\)/2 

If we apply Theorem 15.11 to the stochastic differential equation 



Z t = z o + Ej=i Jo °ij( Z u)dWl in case (a). 

Z\ =Zl- J* (| E U a U Z ")) d u + E U fo °iAZu)dWi in case (6), 

we obtain a transition density To such that we can define the function 
G G C°°([0,1] x R d ) by 



G(t,y):= / r (l-t,y,Og(£)d£,0<t<l 

so that 

[S + iEjU (^ T (^))^rfk) <?(*,!/) = o 

I - Etr (l Ej=! A + I E^> (^(y)) H S^) G ^ v) = 

We define the function F : E — > R by setting 

x), in case (a), 

log(x)), in case (6), 

where logx = (log(xi), log(xd)), and the operator C by 

d 1 d d 2 

C:= m + 2^ Lkl{x) d^dx- l > 

k,l=l 

where L ki (x) = Ej=i ^Ay'O^Vy^)- Now we have that 

CF{t,x) = on [0,1) x £, 



F(t,x) :-- 
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and Ito's formula implies that 



F(t, X t ) = F(0, X ) + J2 [ 4~ F ^ X u)dXt a.s. t G [0, 

k=1 Jo ox k 

From Theorem 15.11 we get that 

F(t,X t )^ f(Xt) inL 2 ast/l 

and 



1). 



fiX,) = F(0, X ) + V / ^-F(u, X u )dX„ 
k=1 Jo ox k 

5.2. Definition. For /, F and X as above we define 



a.s. 



EE 



_d_ 

dxi 



F(u,X u )--^F(^_ 1 ,X tu ))dX^ 



L 2 



for all t = (ti)'i =0 G T n and s G [0, 1). 

5.3. Definition. We define H x f, H* x f : [0, 1) -> [0, oo) by setting 



H x f{t) := ( supE 



Faa(Xt)Lj3p[Xt 



d 2 



dx a dxp 



-F{t,X t 



and 



H* x f{t) := sup H x f(s). 

se[o,t] 



Finally, we define functions Qi : lR d — > K for z = 1, d by 

1, in case (a) 
Xj in case (b). 



Qi(x) :-- 



In this setting we have the following theorem, which refines P, The- 
orem 1]. 

5.4. Theorem. Assume that for all x G E 



QS 



;// Mr a vi X ) 



Qi(x) 



where q + r = s, q, r, s G {0, 1, 2}, 



\bi{x)\ < CiQi(x) and La(x) > -^rQ 2 (x) fori G {l,...,d} and some 
fixed C\ > 0. 

(1) If one has that 

i 



then 



I H := / H x f(t)dt < oo, 



inf sup 4 m (/(X!), r, s) < /or all nGN, 

reT » se[o,i] V ra 



22 



HEIKKI SEPPALA 



where D 1 = D 1 (C 1 , d) > 0. 
(2) // there exists C 2 > and a G (|, 1) such that 

(l-log(l-f))- 



H* x f(t) < C 2 - 

i — i 

then 

inf sup ar(/(X 1 ),r, S ) < 



Do 



se[o,i] 



2a- 1 



/or a/Zt G [0,1), 



/or all nGN, 



where D 2 = D 2 (C\, C 2 , a, c?) > 0. 

Proof of Theorem 15.41 . Hujo showed in the proof of P, Theorem 1, p. 
18] that under the assumptions of Theorem 15.41 we have that 

2 

9 -F(u,X u ) - ^-FitU.X^)] dX k u 



E 



n d 



EE 

1 = 1 k: 

n 

^E 



t "As 



, As 



i=l k=l " 

sup E 



OXk OXk 



& 2 



dx a dxp 



-F(u,X„ 



dudt 



pi 

<cJ2 {ti-t)[H* x f{t)Ydt 
.■ 1 Jtv , 



i=l "H-i 

for any s G [0, 1) and any time net r = (£™)™ =0 , where c 
Hence we can conclude by Theorems 12.31 and 12.51 



c(C u d). 
□ 
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